The fate of chaotic binaries 
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A typical stellar mass black hole with a lighter companion is shown to succumb to a chaotic precession 
of the orbital plane. As a result, the optimal candidates for the direct detection of gravitational 
waves by Earth based interferometers can show irregular modulation of the waveform during the last 
orbits before plunge. The precession and the subsequent modulation of the gravitational radiation 
depends on the mass ratio, eccentricity, and spins. The smaller the mass of the companion, the more 
prominent the effect of the precession. The most important parameters are the spin magnitudes and 
misalignments. If the spins are small and nearly aligned with the orbital angular momentum, then 
there will be no chaotic precession while increasing both the spin magnitudes and misalignments 
increases the erratic precession. A large eccentricity can be induced by large, misaligned spins but 
does not seem to be required for chaos. An irregular precession of the orbital plane will generate 
irregularities in the gravitational wave frequency but may have a lesser effect on the total number 
of cycles observed. 
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Merging black hole binaries are potent sources of grav- 
itational waves and are among the most promising tar- 
gets for direct detection by the future interferometric 
observatories. Black hole mergers, if sufficiently abun- 
dant, are likely to be the most common compact binary 
merger to be detected. If the black holes are rapidly 
spinning, then the orbit can be extremely irregular, even 
chaotic, bearing significant implications for gravitational 
wave searches . An earlier Letter || identified chaos 
in relativistic, spinning binaries in a somewhat abstract 
discussion. In this article, the intention is instead to pro- 
vide a more concrete discussion with less emphasis on 
formal chaos. What is observationally important is visi- 
bly irregular motion. Taking this attitude, a specific as- 
trophysical model is followed through successive stages in 
order to gauge when irregular motion will occur within 
the LIGO/VIRGO bandwidth. Specifically, we investi- 
gate the orbits of a maximally spinning black hole with 
a lighter companion. 

Certain binary star systems are fated to evolve into 
black hole binaries. The orbits of these long lived bi- 
naries have sufficient time to circularize before entering 
the LIGO/VIRGO bandwidth as angular momentum is 
lost to gravitational waves. An archive of circular tem- 
plates is accruing for various binary parameters. Yet the 
merger rates of these evolved binaries are predicted to 
be too low to ensure detection by the first two gener- 
ations of LIGO detectors. A more promising detection 
rate is predicted for dynamical binary black holes; that 
is, binaries formed by the dynamical capture of one black 
hole by another in dense stellar systems Q. The merger 
rate is expected to be about 1.6 x 10~ 7 /yr/Mpc 3 , which 
exceeds neutron star merger rates as well. Dynamically 
formed black hole binaries should have a distribution of 
eccentricities and short orbital periods with masses in 
the range of ~ (5 — 15)M Q M. A binary with masses 



77ii = 15Mq and 777.2 = 5M@ for instance will emit gravi- 
tational waves with a frequency within the optimal LIGO 
bandwidth of / ~ (10 — 10 2 ) Hz for radial separations 
r < 50m where units of total mass m = mi + 7772 are 
used. These provide natural values for the mass and radii 
ranges to investigate. The heavier black hole is taken to 
have maximal spin Si = m\ (spin period P~3x 10~ 4 s 
for a 15Mq black hole). Unlike pulsars, black holes are 
expected to essentially maintain the spin they are born 
with H through most of the inspiral. This canonical 
BH/BH pair can precess chaotically any time the tra- 
jectory transits near the underlying homoclinic orbits of 
Refs. jnj. Homoclinic orbits are purely relativistic, a 
consequence of nonlinearity, and unstable. They have 
the essential features for the onset of chaos |^|,[l0[ when 
the bodies spin. Still, having said this, it is not clear that 
chaos will be confined to this region of phase space. 

A BH/NS binary with typical parameter values of 
777i = 10M Q and 777.2 = 1.4A/q follows trends similar to 
the BH/BH binaries. The explosive evolution of stellar 
progenitors which populate BH/NS pairs delivers large 
kicks to the objects and leads to large spin misalignments 
jnj. It is still unclear whether the population of such 
pairs is too sparse to expect a detection. Since it is only 
the mass ratio which enters the equations, either of these 
cases can be scaled to represent the dynamics of much 
more massive systems which will be visible to LISA [ p"2| . 

We look for chaotic behaviour when energy is con- 
served and the radiation reaction is turned off. For an in- 
dividual orbit, chaos manifests itself as the unpredictable 
precession of the orbital plane. When considering all 
possible orbits, chaos manifests itself as an extreme sen- 
sitivity of the orbital precession to initial conditions so 
two neighboring orbits may live out very different pre- 
cessional histories. The implication is that there is a 
theoretical limit on how well we can predict the orbit 



and therefore the waveform of the emitted gravitational 
radiation Jl],|^]. Dissipation from the emission of gravi- 
tational waves is then included. Irregular motion in the 
dissipative system is understood in terms of the number 
of windings the pair executes in a region of phase space 
which is chaotic for the underlying conservative system. 

The regularity of the orbit will be effected by several 
parameters: the mass ratio (3 = m 2 /mx, the magnitudes 
of the spins, the spin alignment with respect to the orbital 
plane, the eccentricity of the orbit, and the radius of the 
orbit at the time of detection. As is already clear from 
Rcf . Q , motion in the conservative system becomes more 
irregular the larger the angle the spin makes with the 
perpendicular to the orbital plane. The importance of 
three other parameters is evaluated here: (1) the binary 
mass ratio (3 = m, 2 jm,\, (2) the magnitude of the second 
spin S%, and (3) the eccentricity of the orbit. 

The conclusions in brief for the three parameters var- 
ied below are the following: (1) The mass ratio primarily 
effects the cone of precession. The smaller the mass ratio 
P = rri2 /mi, the larger the angle subtended by the orbital 
plane and the larger the modulation of the gravitational 
waves |L3|,|l4j. (2) There can be chaotic motion when the 
massive black hole spins rapidly even if the companion 
has no spin. Still, the larger the magnitude of the sec- 
ond spin (and the misalignment), the more irregular the 
motion. (3) Eccentricity is a consequence of large, mis- 
aligned spins and therefore it is difficult to separate cause 
and effect here. Still, it is clear that eccentricity alone is 
not responsible for chaos. 



might even be inclined to take the extreme resistance of 
the relativistic two-body problem to solution as evidence, 
or at least confirmation, of nonintcgrability. 

The validity of the PN expansion is not questioned 
further and the equations are treated as a self-contained 
dynamical system. In the PN scheme, the orbit evolves 
according to the force equation 



r = a PN + a S o + a S s + a RR 



(2.1) 



in center of mass harmonic coordinates [ jl4fl . The accel- 
eration is due to Post-Newtonian (PN) effects, spin-orbit 
(SO) and spin-spin (SS) coupling and radiative reaction 
(RR). The explicit form of a can be found in the ap- 
pendix. The spins also precess due to the relativistic 
frame dragging and Lens-Thirring effect. The precession 
equations are 
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and 
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s 2 = n 2 x S 2 (2.2) 



L N -S 2 + 3(n-S 2 )n , (2.3) 



L N - Si +3(n- Si)n . (2.4) 



The spins precess with constant magnitude although the 
total spin S = Si +S 2 may not have constant magnitude. 



II. EQUATIONS OF MOTION AND SPIN 
PRECESSION 

The Post- Newtonian (PN) expansion of the relativistic 
two body problem leads to a system of equations describ- 
ing the fate of spinning binaries The PN expansion 
converges slowly to the fully relativistic description [|l6| . 
For this reason, it is a poor approximation at small sepa- 
rations. Despite its shortcomings, the PN expansion does 
give the qualitative features of a relativistic system such 
as nonlinearity, the existence of unstable circular orbits 
p6| , homoclinic orbits JlO|, and spin precession. Since 
these are the ingredients for chaotic dynamics, the qual- 
itative behaviour should persist in a more accurate ap- 
proximation, although the quantitative conclusions are 
subject to change (see for instance the improved tech- 
nique of Ref. Q). 

It is worth emphasizing that approximations can in- 
troduce chaos when the exact system is truly regular. 
One might worry that the error at 2PN order has intro- 
duced chaos which would be removed if we knew the full 
equations of motion without approximation. However, 
the relativistic two-body problem is likely to be more ir- 
regular at higher orders as the nonlinearities of general 
relativity are more accurately represented, not less. One 




orbital plane 
initially 



FIG. 1. A schematic drawing of the inclination angle 
i = arccos (Ln • J) and the angle 9s = arccos (Ln ■ S) . The 
orbital plane traces out a band as in figs. I and I as the 
Newtonian angular momentum precesses about J. 

The orientation of the orbital plane is defined by the 
Newtonian orbital angular momentum 

L N = /z(r x v) (2.5) 

with the reduced mass [i = m\ra 2 lm and the total mass 
m = mi + m 2 . Spin precession, generates a precession of 
the orbital plane (fig. |l|). This can most easily be seen by 



noting that to 2PN-order the total angular momentum J 
is conserved with J = L + S. The orbital angular mo- 
mentum L can be split into two pieces, L = qLn + Lso 
as in eqns. ( A15| ) and (A19). The term Lso is due to 
spin-orbit coupling and a contains Post-Newtonian cor- 
rections. To 2PN order J = and so 

L - -S. (2.6) 

The magnitude of the orbital angular momentum is not 
constant. Further, the precession of the orbital plane 
can be much more complicated than the precession of S: 

^o;Ln) ~ — S — Lso- The orbital plane therefore does 
not just carve out a simple cone as it precesses around the 
direction of J. Instead the plane tilts back and forth as 
it precesses. The motion becomes chaotic as the tilting 
back and forth becomes highly irregular in the Hamil- 
tonian system. When radiation reaction is included, the 
degree of irregularity in the precession of the orbital plane 
can be estimated in terms of how many windings the pair 
spends near the chaotic region of the underlying conser- 
vative system. 

Fig. g shows an example of the trajectory for the cen- 
ter of mass of a BH/NS pair. The motion is plotted in 
three-dimensions to illustrate the precession of the orbital 
plane. 




FIG. 2. A three-dimensional view of a regular orbit with 
/3 — m%l<m\ = 1.4/10, Si = and S2 = 0. The initial angle 
61 = arccos(LN ■ Si) = 45°. The initial conditions for the 
orbit are r/m = 20 and r<j> — 0.209. 

The precession, whether regular or irregular, modu- 
lates the emitted gravitational waves. The response of 
a detector on Earth to an impinging gravitational wave 
can be parametrized as 



where h + and h x are the two gravitational wave polar- 
ization states and F+ and F x are the antenna patterns. 
A radiation coordinate system can be defined such that 
the polarization axes are fixed even in the presence of 
precession [18[ . In such a coordinate system, F + and F x 
are constant. However, the polarization states h + and 
h x depend on the inclination of the orbit and the preces- 
sional frequency in a complicated way [fl4| . For a circular 
orbit, the detector response can be written 



h = Acos(2<f>-5), 



(2.8) 



where the higher harmonics have been ignored for sim- 
plicity. The amplitude A and polarization phase 8 de- 
pend on the binary's location, orientation, and precession 
|l3],[l4j]. For an elliptical orbit, h has terms of the form 
cos($), cos(2$), and 008(3$), at quadrupole order so the 
gravitational wave spectrum shows oscillations at once, 
twice and three times the orbital frequency. 

Precession of the orb ital plane will (1) modulate the 
amplitude in eqn. 



(2-' 



(2) modulate the polarization 
phase 6 and therefore the frequency of the gravitational 
waves and (3) contribute to the overall accumulated 
phase by changing the inspiral lifetime. Any extreme 
sensitivity to initial conditions will most likely have the 
largest effect on the modulation of the amplitude and fre- 
quency of the gravitational waves. The overall accumu- 
lated contribution to the number of cycles in the observed 
waveform will certainly be effected by the general bulk 
precession but may be less sensitive to the irregularity of 
the precession until the very final stages of coalescence. 
The reason for this is that at the radii accessible to the 
interferometers, the irregularity seems to predominately 
effect the orientation of the orbit with a lesser effect on 
the net orbital velocity u = 88 + sin0#. LIGO/ VIRGO 
aim to observe gravitational waves by accurately measur- 
ing the accumulated phase defined as 



* = 



ivdt = I —doj, 



(2.9) 



F x h> 



(2.7) 



which may only have a small correction from the irregu- 
larity of the precession. Irregular motion will effect the 
phasing, that is the gravitational wave frequency, and the 
amplitude of the wave. The number of cycles in the accu- 
mulated phase may therefore be determined by the aver- 
age bulk behaviour of the precession. Although this con- 
jecture requires further scrutiny, we can use the approxi- 
mations of Refs. Jl^,0 to provide a rough figure for the 
number of wave cycles. In Ref. Jl3|, circular Newtonian 
orbits were studied. The effect of spin precession was 
isolated without including spin-orbit acceleration terms 
in the equations of motion in Ref. |l3| . This separa- 
tion of the spin precession equations and the equations 
of motion removed any possibility of chaotic coupling. 
However, the gross features can be fairly represented by 
this approach. They estimate that the change in the to- 
tal number of cycles amounts to about twice the total 
number of precessions. For the binary black holes of size 



15M Q + 5Mq, there can be ~ 10 — 15 precessions in the 
observable band and so there can be ~ 20 — 30 addi- 
tional cycles (depending on spin orientation, magnitude 
and eccentricity of the orbit). In Ref. fl4[ , the contribu- 
tion to the total number of orbits over the entire LIGO 
bandwidth was also estimated to be ~ O(20) for black 
hole binaries with total mass ~ 20Mq. The number of 
additional cycles could quadruple for NS/NS masses |lj] 
It has been argued that matched filtering will be a 
poor means by which to observe BH/BH orbits in the 
chaotic regime PJl9| and that other cruder methods will 
be employed. We contribute to this debate only by in- 
dicating when irregularity will influence dctectability. It 
has also been emphasized in Ref. [jlij, that the PN ex- 
pansion is not sufficient to accurately model templates 
for r < 15m. Nonetheless, the higher order contributions 
will incorporate stronger relativistic effects and therefore 
more nonlinearity which should only lead to more irreg- 
ular motion. Since we are trying to provide a physical 
picture of the general trends and dependences, we con- 
tinue to use the PN expansion. 



III. CHAOS 

Chaos in relativistic systems is notoriously difficult to 
quantify |^0|,|l| . While formal definitions of chaos are of 
little interest to data analysts when irregularity regard- 
less poses a hindrance, the tools of chaos are nonetheless 
critical to survey the system for endemic irregularity. In 
Ref. p|, chaos in the spinning black hole problem was 
identified through the method of fractal basin boundaries 
[^[U^] . In this section irregularity of an individual orbit 
is discussed and the power of the fractal basin boundary 
technique is utilized. It is worth emphasizing that the 
fractal basin boundary method allows one to scan huge 
numbers of orbits and therefore provides an invaluable 
tool to survey phase space for chaos. 

The equations of motion and spin precession equations 
are treated as any other dynamical system. In the first 
instance radiation reaction is omitted and chaos is han- 
dled in the conservative system. Dissipation is treated 
in § VII. The 2PN equations of motion can be derived 
from a Lagrangian [^2|,^3) and therefore can be derived 
from a Hamiltonian. The coordinates Si(i) and S2(t) 
can be treated as external time-dependent perturbations 
to the integrable system with the eq uatio ns of motion 
supplemented by the precession eqns. The system 

can therefore be treated as a time-dependent Hamilto- 
nian system H(r, v, Si (i) , S2 (t)) . 

Chaos is well defined for Hamiltonian systems. Chaos 
is synonymous with nonintegrability. Regularity is syn- 
onymous with integrability. In a Hamiltonian system 
with N degrees of freedom q and N conjugate momenta 
p, integrability prevails when there are N independent 
constants of the motion. The constants of motion must 
also be in involution; that is, the Poisson bracket of any 



constant with the others vanishes: [Cj, Cj] — 0. The 2N- 
dimensional phase space can then be reduced to motion 
on an TV-dimensional torus. This is most easily repre- 
sented with a canonical change of coordinates to action 
angle variables (I, 0) such that each of the N momenta 
/ is set equal to one of the N constants of motion C . The 
motion can then be made periodic in the coordinate vari- 
able 8. For N = 1, degree of freedom, the motion lies on 
a circle. For TV = 2, the motion lies on a torus and for 
arbitrary N, the motion lies on an TV-dimensional torus 
^4j,^5|. In any other set of canonical variables besides 
action angle variables, the mark of integrability is that 
the motion is confined to a smooth closed curve in (p, q) 
and does not diffuse off that line. If the motion in (r, f) 
has diffused off of a smooth line it is not restricted to a 
torus and the motion must be nonintegrable. 





FIG. 3. A regular orbit with [3 — m^/mi — 1/3, Si = m\ 
and S2 = 0. The initial angle 9\ — arccos(LN-Si) = 45°. The 
initial conditions for the orbit are r/m = 20 and r(j) = 0.209. 
Top: three-dimensional view of the orbit. Middle: a projec- 
tion of the orbit onto the (r, r) plane. 

For coalescing binaries, there are N = 3 degrees of 
freedom, (r,6, 4>). When there are no spins, the phase 
space is 27V = 6-dimensional. The energy and angular 



momenta provide enough constants of the motion to re- 
strict the trajectories to tori and there is no chaos to 
2PN order Beyond 2PN order on the other hand, 

the two-body problem is likely to be chaotic even in the 
absence of spin. 

When the bodies spin, the dynamics can be reduced 
to a time-independent Hamiltonian system by taking a 
Poincare surface of section. This method of identifying 
chaos is less than ideal here. For the record, we will take 
the time to explain the method and its shortcomings. 
The Poincare surface of section involves plotting a point 
in (r, f) each time the orbit crosses a surface on which 
all of the other coordinates are fixed. If the collection of 
points defines a smooth curve, then the motion is con- 
fined to a torus and is integrable. If the points speckle 
the plane, then the motion has diffused off of a torus 
and is nonintegrable. The Poincare reduction of phase 
space is most easily demonstrated for the case of only 
one spinning body so that S = Si. Including spin, the 
phase space is 2N + 3 = 9 dimensional. The 4 constants 
of motion, E, J, restricts motion to a 5-dimensional sur- 
face. If we treat the Hamiltonian as periodic in the three 
coordinates S(t), we can plot a point in (r,r) each time 
the orbit crosses S(t) = S(0) so that there are only two 
free coordinates remaining. In this way we can look for 
the destruction of tori and test for nonintegrability. 

A regular BH/BH binary with mass ratio (3 = 
m,ijm\ = 1/3 is shown in fig. |3| The more massive black 
hole has spin S\ = rri\ with arccos(LN • Si) = 45° ini- 
tially. The second spin is set to zero for this figure. The 
top panel in fig. ^ shows a three-dimensional view of the 
orbit. The middle panel shows a projection of the or- 
bit in the (r, f) plane. This orbit is very nearly circular. 
Note that if the orbit is exactly circular with only one 
spin then there can be no chaos since the dimensionality 
of the phase space reduces to two , which is not enough 
to support chaos. This orbit is regular as indicated by 
the absence of spreading in the projection onto the (r, f) 
plane. 

By contrast, chaotic orbits are shown in figs. |J and ||. 
Chaos is first identified by the method of fractal basin 
boundaries. To build the basin boundary, 40,000 black 
hole binaries are evolved varying only the initial velocities 
(f, r<f>). The initial condition is then color coded accord- 
ing to outcome: black if the pair merge and white if the 
pair execute at least 50 or more orbits. (Increasing the 
required number of orbits will shrink the white basins but 
will not eliminate structure at the boundaries.) A frac- 
tal boundary signifies extreme sensitivity to initial con- 
ditions and a mixing of orbits - i.e., chaos. The power of 
the fractal basin boundary method as a survey of large 
collections of orbits is clear. Orbits near the boundary 
will be chaotic. 




FIG. 4. A fractal basin boundary scan in (f, nj>) 
varying over initial values in < r < 0.035 and 
0.425 < r<j> < 0.443125. All 40,000 orbits have a maximally 
spinning black hole and a lighter companion with no spin 
(/3 = mijm\ = 1/3, Si = m\ and S2 =0). The initial angle 
is 6\ — arccos(I/N ■ Si) = 95° and initial separation r/m — 5. 
200 x 200 orbits shown. 
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FIG. 5. Top: A three-dimensional plot of an irregular 
orbit near the basin boundary of fig. ^ with initial velocities 
r = 0.00105 and r<f) = 0.43075. Bottom: A projection of the 
orbit onto the (r,r). 



An irregular orbit selected near the fractal boundary 



is shown the top panel of fig. The bottom panel of 
fig. |^ is a detail of the projected motion in the (r, f) 
plane which shows some threading of the orbit. A full 
Poincare section taken as described above confirms that 
this is not an illusion from the projection. However, the 
numerical burden is extreme and impractical in nearly all 
cases of interest. For this reason, we continue to rely on 
the fractal basin boundaries and only take the projection 
onto the (r, r) plane as a crude guide and not as proof of 
chaos. 

In theory, the surface of section technique can also be 
implemented in the case of two spinning bodies since both 
are periodic under the precession. In practice this can be 
difficult since one might have to wait a very long time 
before Si(i) = Si(0) at the same time as S 2 (i) = S 2 (0). 
Alternatively, we can cheat and simply look at the pro- 
jection of (r, r) in the full phase space without taking a 
Poincare section. This is only used as a crude survey tool. 
If the projected motion lies on a simple closed curve, the 
motion is decidedly integrable. If the motion lies on a 
threaded orbit (such as that of fig. ||) , then this indicates 
the motion might have diffused off of a torus. 

Another major shortcoming of the Poincare surface of 
section method in this setting is that the 2PN constants 
of motion are only approximately conserved. Therefore 
even if one is cautious, the spreading may be a result of 
the approximation and not a true mark of the destruction 
of tori. Since this projection is ambiguous, it is only used 
as a rough guide. For a firm identification of chaos we 
rely on the unambiguous fractal basin boundaries. 

There are other outcomes that could be used to color 
code a plane in phase space and study the basin bound- 
aries. A set of outcomes based on the number of windings 
a pair execute would be relevant to the search for gravita- 
tional waves. In fig. ^|, basin boundaries for the same set 
of initial conditions are compared. The binaries all have 
the same initial condition except for the initial angular 
displacement of the spins. The top panel was originally 
published in For that top panel, stable/merger out- 
comes are used. The initial condition in the (8\, 82) plane 
is color coded white if the pair executes at least 50 orbits 
and black if the pair merge in under 50 orbits. The lower 
panel uses a winding criterion: the initial condition in 
the (61, 82) plane is color coded white if the pair execute 
more than 50 orbits before merger, black if they merge 
after executing more than 36 orbits (but less than 50), 
dark grey if the pair merge having executed 36 orbits, and 



*A fruitful analytic approach may be to treat the motion 
as an instance of Arnol'd diffusion H|. Even if r is small 
and S changes slowly, their presence induces a coupling be- 
tween 8 and <fi which can lead to chaotic resonances and hence 
stochastic behaviour. The future direction of this work is to 
interpret the chaotic motion in terms of this slow modulation 
diffusion @. 



light grey if the pair merge having executed less than 36 
orbits. The mixed basin boundaries show that there is 
some unpredictability in the number of orbits executed 
in the conservative system. 




FIG. 6. All orbits have m-z/mi — 1/3, r — 5m, and 
r<j) = 0.45. The initial angles are varied in the (81,02) 
plane over the ranges 140.375° < 6>i < 157.5635° and 
45.8367° < 8 2 < 63.025411°. Top: Figure from The frac- 
tal basin boundaries built up using the stable/merger crite- 
ria. Bottom: The fractal basin boundaries built up using a 
critical winding criterion. A resolution of 200 x 200 orbits. 

We discuss dissipation in jjvT|. We point out here that 
when dissipation is included, there is less fractal struc- 
ture when the critical windings are used as outcomes. 
This indicates that with dissipation included the number 
of windings executed by a binary might be predictable. 
However there is still structure in the basin boundaries 
when other outcomes are used indicating that not all fea- 
tures of the wave form will evolve predictably. 



IV. THE BINARY MASS RATIO 

The binary mass ratio (3 = 777,2/7771 primarily effects the 
cone of precession. The lighter the relative mass of the 
companion, the larger the band in which the orbital plane 
will precess and the l arg er the corresponding modulation 
of the waveform fl^fbf . The orbital plane will precess 
around J with an angle of inclination 



cos i = Ln • J = 



V 7 ! + (s/L N y 



(S/L N )cos6 s 



2{S/L N )cos6 s 

(4.1) 



where cos6>s = Ln • S as in fig. [I] and the total angular 
momentum has been used to lowest order, J ~ Ln + S. 
The angle subtended changes as the angle between Ln 
and S changes. This leads to the additional tilting back 
and forth on top of the simple precession. The ratio 
S/Ln can be estimated by taking Ln ~ fi(mr) 1 / 2 for a 
nearly circular orbit. Consider the extremes when the 
more massive object spins with S\ = m\ to get an upper 
range and the opposite regime when the lighter star spins 
with 5*2 = 7712 to get a lower range: 



777 2 



S 



(m/r) 1/2 < — < 
7771 / Ln 



777l 
7772 



(m/r) 



1/2 



(4.2) 



For very small S/Ln, the precession cone angle is so tight 
that modulations of the gravitational wave signal will not 
be s ubstantial. Notice from the left-hand side of eqn. 
(4.2) that if only the lighter object spins then S/Ln is 
small for all r. This may explain why Ref. jjj found 
chaos for a test particle around a Schwarzschild black hole 
only if the light companion had a spin 5*2 > 0.64m2mi 
for mi 3> 7772, which exceeds the maximal spin of S2 = 
m|. The 2PN expansion to the two-body problem allows 
both black holes to spin and precess. When the more 
massive object spins maximally, then S/Ln is large for 
mi 3> 777 2 . For a 10 : 1 mass ratio, S/Ln ~ 1 for r ~ 
IOOto. As a result, for large mass ratios, the orbital plane 
subtends a larger angle at a given radius and modulation 
will be correspondingly larger [ p^|Jl4| . Consider fig. || 
versus fig. In both figures, the more massive BH spins 
maximally [Si — to 2 ) with a spin displacement of Ln • 
S = cos(45°) = l/y2. In both figures, the companion 
has no spin. The difference between the two figures is 
the mass ratio. In fig. |^, the mass ratio is (3 — 1.4/10 
and the angle subtended at r — 20 is i ~ 28°. In fig. 

the mass ratio is (3 = 1/3 and the angle subtended at 
r = 20 is i ~ 18°. The band subtended by the precessing 
plane is larger for the smaller mass ratio, although it is 
substantial in both cases. 



We can invert eqn. (4.1) using the right-hand-side of 
eqn. (|4.2|) to write the radius at a given inclination, 



, r v 1/2 
V777/ 



777l 
7772 



COS 
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with cos 81 = Ln ■ Si. We could take this as an indicator 
for the radius at which precession becomes important. 
Letting i ~ 15° and 9\ ~ 45°, then 
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(4.3) 



For circular motion, L ~ /i(mr) 1 / 2 and m(j) ~ (m/r) 3 ^ 2 
and the frequency of the emitted wave, / ~ cft/n is 
roughly 



3fl3 



/ - 9 x 10* ft 



Me 
m 



Hz. 



(4.4) 



For a BH/BH binary with (3 ~ 1/3 and 777 = 20M Q , then 
/ ~ O(10) Hz when the spin precession angle opens to 
7 ~ 15° and the effects on the gravitational wave should 
be noticeable. For a BH/NS binary with (3 ~ 1.4/10 and 
777, = 11.4, then / ~ 0(\) Hz when the spin precession 
angle opens to i ~ 15°. To emphasize, precession will be 
important for these and larger frequencies as the pairs 
sweep through the interferometer's bandwidth. 
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FIG. 7. The waveform, Newtonian angular momentum and 
spin for the orbit of fig. H. Top: The waveform h+. Bottom: 
three-dimensional view of the precession of Si. 



While we can conclude that smaller mass ratio means 
a wider precessional angle, the dependence of the reg- 
ularity of the motion on (3 within that band is not yet 
clear. Within the bulk precession, the orbital plane tilts 
back and forth, sometimes regularly and sometimes er- 
ratically. Orbits can subtend roughly the same cone but 
some occupy the band more regularly than others. A hint 
of the effect of (3 on regularity comes from a comparison 
of frequencies. An instantaneous orbital frequency can 
be defined roughly as 



Ln 

2 



(4.5) 



for comparison with the instantaneous spin precession 
frequency of the larger object (neglecting spin-spin cou- 
pling just for the crude estimate) 



so that 



3m : 
2mi 



0(2 + 3/3/2) Vm 
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(4.7) 



Therefore uj/Q,\ is always large. The pair executes many 
orbital windings per spin precession. This hints that the 
chaotic behaviour may be related to modulation diffusion 
where a slowly varying parameter, in this case Si, facili- 
tates chaotic resonances between coordinates, in this case 
6 and $ @. 

However, compared to the instantaneous precession 
frequency of the spin of the lighter object 



n 2 



09 + if 
(20 + 3/2) Vm 



(-) 

Vm/ 



(4.8) 



This ratio is not as large for small (3 so the lighter com- 
panion always precesses much more than the heavy black 
hole. As argued in §|v|, the spin of the companion en- 
courages irregularity and the fact that 17 2 /Sli < 1 may 
account for some of this effect. 

In general, the mass ratio f3 = 7712/7711 seems to effect 
the bulk shape of the prec ession and gravitational wave 
modulation. By eqns. (4/7) and (O), [3 may determine 
the radii at which chaotic resonances will occur. Still, it 
is unclear how much the mass ratio impacts the details 
of the motion. 

For comparison with later cases, the waveform emitted 
by the BH/BH binary of fig. || is shown in the top panel of 
fig. f^. Even though the precession is fairly regular, it does 
modulate the wave amplitude and frequency. The mod- 
ulation of the waveform has been minimized by placing 
the binary directly above the detector. Tail-effects are 
neglected throughout. A three dimensional view of the 
precession of the spin Si is shown in the lower panel of 



V. SECOND SPIN 

The effect of spinning up the lighter companion can 
be studied by starting with the orbits of fig. |J. Even 
a small second spin will cause further diffusion in phase 
space. When the second spin is maximal, the fraying 
of the projection in (r, r) increases. Fig. || shows an ir- 
regular three-dimensional orbit, the projected motion in 
(r, f), and the waveform when both objects spin maxi- 
mally. Fig. U shows the precession of both spin vectors. 
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FIG. 8. An irregular orbit with /3 — m%lmx = 1/3, 
Si = m\. The initial conditions for the orbit are r/m — 6 
and f — 0.025 r<j> = 0.365. Top: three-dimensional view of 
the orbit. Middle: a projection of the orbit onto the (r, r) 
plane. Bottom: the waveform. 



seems to occur when angular momentum sloshes between 
spins and the orbit. It is difficult to separate cause and 
effect. Still, it is obvious that eccentricity alone is not 
the culprit. 



The fraying of the orbit in (r, f) indicates that the mo- 
tion might not be confined to a torus. As discussed in 
§111, the projection is a hint of nointegrability; that is, 
of chaos. However, as already mentioned, the projec- 
tion alone is not enough to conclude there is chaos. A 
full basin boundary analysis shows that this orbit occurs 
near a fractal boundary giving every indication the orbit 
is chaotic. Notice in fig. [9] that the spin of the lighter 
star precesses more tha n th e spi n of the heavier object 
as expected from eqns. (|4.7|)-(^8|). 





FIG. 9. The Spins for the orbit of fig. 
three-dimensional view of the precession of Si. 
three-dimensional view of the precession of S2. 



Top: 
Bottom: 





VI. ECCENTRICITY 

Binary black holes formed in dense stellar regions are 
thought to have a roughly thermal distribution in eccen- 
tricity with a slight enhancement of high e at the time of 
formation While many of these may still have time 
to circularize before merger, it is worth investigating the 
role of eccentricity on the regularity of the orbit. Large 
misaligned spins necessarily induce eccentricity. Chaos 
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FIG. 10. A regular but eccentric orbit with 
f3 = 7712/mi = 1/3, Si — rrix and S2 = 0. The initial an- 
gle 9i = arccos(LN • Si) = 45°. The initial conditions for the 
orbit are r/m = 20 and r<j> — 0.2. Top: three-dimensional 
view of the orbit. Middle: a projection of the orbit onto the 
(r, f) plane. Bottom: the waveform. 



Consider fig. [lC] where only one of the holes is spin- 
ning. The top panel shows the three-dimensional orbit 
which does have irregular features. However, this is de- 
ceptive. It is well known that Keplerian orbits are closed 
ellipses while relativistic, elliptical orbits precess within 
the orbital plane. The entire plane then precesses due 
to the spins. What is being witnessed in the top panel 
of fig. [lO] is this double precession. The middle panel 
shows complete regularity of the orbit in (r, r). Motion 
in this coordinate is confined to a torus. The waveform 
shown in the bottom panel shows several expected fea- 
tures. Since the orbit is elliptical, the gravitational waves 
oscillate at once, twice and three times the orbital fre- 
quency which changes the spectrum from that for a cir- 
cular orbit. The double precession then modulates the 
amplitude and phase on top of these oscillations. Even 
though this motion is regular (for S2 = 0), in the sense of 
being predictable, the modulation of the waveform from 
the double precession must certainly impact observations 
gained through the method of matched filtering. 

Eccentric orbits do show chaotic precession when the 
companion spins rapidly as well. (Note that when both 
stars spin, there are no circular orbits [B.) 

VII. DISSIPATION 

The emission of gravitational waves dissipates energy 
and angular momentum. The spins are not strongly ef- 
fected by the loss of gravitational waves which carry away 
orbital angular momentum and circularize the orbit. The 
irregularity of a dissipative binary can be evaluated in 
terms of how many orbits the pair executes as it passes 
through the successive regions in the conservative system 
p6[ . The most efficient way to do this is again using the 
fractal basin boundaries. As pointed out in Ref. p6| , 
when dissipation is included the boundaries will never 
be truly fractal. Like a snowflake, the self-similar struc- 
ture will cutoff at some physical limit. However, if a 
color coded basin boundary shows substantial structure, 
it is fair to say, the participating orbits are influenced 
by irregularity. The smaller the scale at which the cut- 
off occurs, the more irregular the history of that set of 
orbits. 

Fig. [ll] illustrates the dissipative inspiral of 90,000 
maximally spinning black hole pairs with (3 — 1/3. The 
orbits all coalesce due to energy lost to gravitational 
waves. The initial condition in (#1,6*2) is color coded 
black if the pair merges from above the z-axis and white 
if the pair merges from below the z-axis The dissi- 
pative system does not provide any ideal outcomes, but 
this criterion should not falsely introduce structure. The 
top panel of fig. [H] shows some structure for maximal 
spins while the lower panel of fig. |ll| shows substantially 
less structure for spins 1/10 of maximal. Incidentally, 
even more structure is seen for BH/NS binaries with 
0=1.4/10 §. 




FIG. 11. Coalescing orbits with radiative reaction for 
/9=l/3 with r/m = 25. The slice through phase spans 
scans over initial orientation of the spins: —180° < 9i < 180° 
and —180° < 82 < 180°. Upper Panel: Spins are max- 
imal Si = mf. Lower Panel: Spins are 10 times smaller, 
Si = 0.1 rrii . 300 x 300 orbits shown in each panel. 

The orbits begin fairly regularly although the preces- 
sion modulates the waveform. Given unlimited access to 
theoretical templates, matched filtering could in princi- 
ple glean a confident signal, at least until the pair drew 
near the unstable orbits. The orbit becomes more irreg- 
ular as the separation closes and the pair passes near the 
underlying conservative trajectory in fig. || just before 
plunge. Since matched filtering relies on the template 
remaining in phase with the data for many cycles, one 
might hope that a disruption of the last few cycles would 
not be serious. Many authors have already argued for the 
use of other methods at these close separations (see [HJ 
and the references therein). However, others have argued 
p7f that these last few orbits are heavily weighted and 
therefore critical for a successful detection. More impor- 
tant than the visually obvious amplitude modulation the 
precession modulates the gravitational wave frequency. 



A frequency space analysis is still required to determine 
if the modulation inhibits detectability of such irregular 
waves. 



VIII. SUMMARY 

It is reasonable to conclude for a typical stellar mass 
BH/BH system that if the primary black hole spins 
rapidly, the orbital plane will precess unpredictably for 
the final orbits before merger. The dependence of the 
precessional motion was tested as a function of three pa- 
rameters: 

• /3 = mijm\: As was already known, the smaller the 
mass ratio, the thicker the band occupied by the process- 
ing orbit ]l3| , |l^ |. A ratio of (3 = 1/3 is small enough for 
precession to modulate observable gravitational waves. 
The mass ratio may also determine the radius at which 
chaotic resonances occur by effecting the relative orbital 
and spin frequencies. 

• spin magnitudes: The transition to chaos occurs as 
the spin magnitudes and misalignments are increased. 
There can be chaotic motion even if only one body spins. 
If it is a light companion spinning, then the magnitude 
of the spin has to be larger than maximal JjJ. If it is 
the heavier object, then chaos can occur for rapid but 
physically allowed black hole spins. This is true for both 
BH/BH and BH/NS pairs. For a given eccentricity and 
radial separation, a transition to chaos can occur as the 
spin of the companion is increased. 

• eccentricity: Large misaligned spins cause eccentric- 
ity which may therefore be a feature common to chaotic 
orbits. However, large eccentricity alone certainly does 
not cause chaos. If chaos is indeed occurring near the 
unstable homoclinic orbits of Ref. [Q, then we might 
expect a fair spread in eccentricities for chaotic orbits 
since the homoclinic orbits occur with a broad range of 
eccentricities. Importantly, even for a completely regular 
orbit the usual relativistic precession of an elliptic orbit 
in the plane becomes superposed on the precession of the 
orbital plane. This combination, though not an indica- 
tion of chaos, does lead to a complicated modulation of 
the gravitational wave signal. 

The most pressing question remains, will irregular 
orbits of coalescing binaries hinder observations? A 
detailed study of the modulated gravitational wave is 
needed. A conjecture is that an erratic precession of the 
orbital plane will result in an erratic frequency or phasing 
of the gravitational wave but will not alter the total num- 
ber of cycles much. If this conjecture is fair, observations 
in the chaotic regime may be salvaged with a modification 
of the matched filtering method. As direct detections be- 
come more acute, we can aspire to watch the transition 
from regular to irregular motion. We might then witness 
the signature nonlinearity of general relativity determine 
the fate of chaotic binaries. 
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APPENDIX A: 2PN EQUATIONS OF MOTION 



In the notation of Ref. 14 1 , the center of mass equations of motion in harmonic coordinates are 



r = a PN + a S o + a S s + a RR 



(Al) 



with r = rf . The right hand side is the sum of the contributions to the relative acceleration from the Post-Newtonian 
(PN) expansion, the spin-orbit (SO) and spin-spin (SS) coupling and from the radiative reaction (RR). The explicit 
terms are quoted from Ref. The following notation is used: v = dr/dt, n = r/r, [i = mim^/m, r\ = /i/m, 

8m = mi — ?Ti2, S = Si + S2, and A = m(S?/m2 — Si/mi). The apN = sn + aiPN + a*2PN with 
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The radiative reaction term is due to terms to 5/2PN order and can be expressed as 
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The spin-orbit acceleration is 
1 



aso 



-> 0111 —* 

6n[(nx v)-(2S + A) 
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and the spin-spin acceleration is 



ass = t ^ fi(Si • S a ) + S^n • S 2 ) + S 2 (n • Si) - 5n(n • S a )(n • S 2 ) 

fir- 



(A7) 



where Epm — -Ejv + ^ipjv + E 2 pn and 



1. Constants of the Motion 

E = E PN + E S o + Ess, 



. 1 9 m 
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The total angular momentum is given by 



where 



with Lpn = Ln + Lipn + L2PN and 



J = L + S, 
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